Structure formation in electromagnetically driven granular media 
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We report structure formation in submonolayers of magnetic microparticles subjected to periodic 
electrostatic and magnetic excitations. Depending on the excitation parameters, we observe the 
formation of a rich variety of structures: clusters, rings, chains, and networks. The growth dynamics 
and shapes of the structures are strongly dependent on the amplitude and frequency of the external 
magnetic field. We find that for pure ac magnetic driving at low densities of particles, the low- 
frequency magnetic excitation favors clusters while high frequency excitation favors chains and 
net-like structures. An abrupt phase transition from chains to a network phase was observed for a 
high density of particles. 

PACS numbers: 45.70.Mg 



Large assemblies of macroscopic particles subject to 
external driving, such as vibration, shear or rotation, ex- 
hibit fascinating collective behaviors. Their dynamics is 
poorly understood, especially when inter-particle interac- 
tions are strongly dissipative 0. Furthermore, additional 
complications arise when the grain size goes below 0.1 
mm, and non-trivial interactions due to charging or mag- 
netization come into play. When small particles acquire 
an electric charge or magnetic momentum, the dynamics 
are governed by the interplay between long-range electro- 
magnetic and short-range contact forces. On the other 
hand, precise electromagnetic excitation can be used to 
control the morphology of the resulting pattern. 

Recently, a number of experimental studies was per- 
formed with vibrofluidized magnetic particles H0. Sev- 
eral interesting phase transitions were reported, in par- 
ticular the formation of dense two-dimensional clusters 
and loose quasi-onedimensional chains and rings. Be- 
sides direct interest to the physics of granular media, 
these studies may provide insight into the fundamental 
problem of dipolar hard sphere fluids where the nature of 
solid/liquid transitions is still debated Alternatively, 
vibrofluidized magnetic particles can be considered as an 
extremely simplified model of a ferrofluid, where similar 
experiments are technically difficult to perform. Previ- 
ous studies were limited to very small number of particles 
(about 1,000) due to the intrinsic limitation of the me- 
chanical vibrofluidization technique. 

We developed a technique to electrostatically drive 
fine conducting powders 0, 0|- Our approach offers 
unique new opportunities compared to traditional vibra- 
tion techniques. It enables one to deal with extremely 
fine powders which are not easily controlled by mechan- 
ical methods. In addition electrostatic driving allows to 
control the ratio between long-range electric forces and 
short-range collisions by changing the amplitude and the 
frequency of the applied electric field. 

In this Letter, we extend our studies towards magnetic 
microparticles. In addition to electrostatic excitation, 
we explore the driving of magnetic particles using an ac 



magnetic field. Our studies reveal a rich diversity of be- 
haviors, including formation of dense immobile clusters, 
as well as quasi-onedimensional chains and rings. We 
provide strong evidence of a first order phase transition 
from finite length chains of particles to infinite networks 
as the driving parameters are varied. 

The experimental setup is shown in Fig. fj^ The de- 
sign of the electrostatically driven cell is similar to that 
reported earlier Conductive particles are placed be- 
tween two horizontal transparent conducting glass plates 
(12x12 cm with the spacing of 1.5 mm). An external 
magnetic field is provided by a set of large magnetic coils 
(30 cm in diameter) placed around the cell. To excite the 
granular media contained in the cell, a voltage of - 2 kV 
with a frequency of - 150 Hz is applied to the plates. 
The magnetic interaction is controlled by supplying a dc 
electric current to one of the coils to create a constant 
magnetic field in the range of 0-80 Oe. In addition, pure 
"magnetic" driving of the particles is carried out by ap- 
plying an ac current to a second coil with small or no 
electric field applied to the cell. The amplitude of the 
ac magnetic field can be varied from 0-15 Oe with the 
frequency from to 300 Hz. The granular media used in 
our experiment consisted of spherical Nickel microparti- 
cles with average size of about 90 (im (Alfa Aesar Com- 
pany) . Magnetic moment per particle at the 80 Oe field is 
1 • 10~ 5 emu; saturated magnetic moment is 2 • 10 -4 emu 
per particle, saturation field is about 4kG. The number of 
particles was varied in the range 120, 000 — 300, 000. Ex- 
periments were also performed with 40 fi sized particles, 
but no qualitative difference was found. The experiments 
were carried out in air and in non-polar low- viscosity liq- 
uid (toluene). Real time images were acquired with a 
high-speed- digital camera. 

Electrostatic driving works as follows. Conducting par- 
ticles acquire an electric charge when they are in contact 
with the bottom conducting plate of the cell. When the 
magnitude of the electric field in the cell exceeds some 
critical value, the upward electric force acting on the 
charged particle overcomes the gravitational force and 
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the particles are driven towards the upper plate. Upon 
contact with the upper plate, the particles recharge and 
fall back to the bottom plate. This process repeats in 
a cyclical fashion. The elevation of the particles off the 
bottom plate can be adjusted by the frequency of the ap- 
plied ac electric field. Pure "magnetic" driving works in 
the following way. Magnetic particles with moment M 
subjected to an external magnetic field H experience a 
torque — [M x H\ forcing the magnetic moment of parti- 
cles to be aligned with the applied magnetic field. There 
are, however, two ways in which the magnetic moment 
of the particle can adjust its direction: (i) the magnetic 
moment can rotate inside the particle against the inter- 
nal magnetic anisotropy field and (ii)the whole particle 
can rotate in order to keep the moment aligned with the 
applied magnetic field. In the latter case, the magnetic 
particle needs to overcome the resistance from the fric- 
tion and adhesion forces between the particle and sur- 
face of the plate. Thus, as the magnitude of the external 
magnetic field exceeds some critical value, each parti- 
cle begins to rotate by keeping its moment aligned with 
the field and moves being driven by magnetic drag force 
F rn = MVH| om i where H; oca / designates the local mag- 
netic field coming from dipolar fields of the neighboring 
particles and the external magnetic field. The effective 
time the particle can move may be controlled by the ad- 
justing the frequency of an ac magnetic field. 

As reported earlier, [^, 0] the particles remain immo- 
bile on the bottom plate of the cell if the applied electric 
field is less than some critical value E\. When the field 
exceeds a second threshold value, E 2 > Ex, the system 
undergoes a transformation to a gas-like state. When de- 
creasing the electric field E below E 2 < E < E 2 ), 
small densely packed clusters begin to nucleate. Clus- 
tering is promoted through electrostatic screening (two 
particles in contact acquire a smaller charge than two 
well-separated particles). In the case of magnetic parti- 
cles, an additional complication arises from the dipole- 



FIG. 2: Phase diagram in the magnetic field, voltage plane 
for 90/iin Nickel particles in the air filled cell. Surface fraction 
of particles is 9% (w 200,000 particles) 



dipole magnetic interactions and the dynamics of cluster 
formation and critical electric fields E\ , E 2 are expected 
to change. The phase diagram delineating the primary 
experimental regimes as a function of applied voltage and 
external magnetic field for 200,000 nickel particles with 
average size of 90 ^m (about 9% monolayer coverage) 
is shown in Fig. [3] The frequency / and amplitude V 
of the ac electric field applied to the cell were chosen 
to provide quasi-two dimensional motion of the particles 
(V<1000V, / = 100Hz). The cell was subjected to an 
external dc magnetic field from to 80 Oe. 

The behavior of the magnetic particles in low exter- 
nal magnetic fields (below 20 Oe) is rather similar to 
that reported for nonmagnetic conductive particles. The 
isolated particles are immobile until the electric field ex- 
ceeds a critical value E\. At electric field values above 
E 2 , the granular medium transforms into a uniform gas- 
like phase. Cluster formation is observed in the interval, 
Ei < E < E 2 , in agreement with earlier studies on non- 
magnetic granular systems 0,0. However, the cluster- 
ing dynamics of the magnetic particles is very sensitive 
to the external magnetic field. Clusters formed in the 
absence of a magnetic field (a) and at a magnetic field 
of 10 Oe (b) are shown in Fig. [21 The cluster grown in 
an external magnetic field looks less compact and suffi- 
ciently less smooth in shape. Moreover, cluster formation 
in an external magnetic field is characterized by arrested 
coarsening due to depletion of the gas phase. 

We analyzed the growth rates to characterize the time 
evolution of the clusters at different magnetic fields. We 
find that the change of the cluster size for both zero and 
nonzero external magnetic fields can be well described by 



a simple equation: 

S(t) = Soc(l-exp(-(t-t )/r)) (1) 

Here Soo is the cluster size at t — > oo; to is the time when 
the cluster nucleates; r is the growth time. The quan- 
tity (1 — S(t)/Sao) for the clusters formed at zero field 
and at 10 Oe is plotted in Fig. Solid lines in Fig. [3] 
represent the functional dependence of Eq. Q with the 
corresponding characteristic growth times. The larger 
growth time r in the presence of an external magnetic 
field could be explained as follows: since the particles 
are large enough to contain multi-domains, the higher 
external magnetic field induces larger magnetic moments 
in the particles (due to the growth of magnetic domains 
within the particle with orientation along the external 
field) leading to amplification of the dipole-dipole inter- 
actions. Consequently, the interaction between particles 
becomes highly anisotropic and increases the character- 
istic growth time since (a) aggregation proceeds predom- 
inantly through coalescence of chain segments which are 
less mobile; (b) some particles are repelled from the clus- 
ter due to unfavorable magnetic moment orientation. For 
elevated magnetic fields (20-45 Oe) the formation of im- 
mobile vertical chains consisting of two to four particles 
were observed when the amplitude of the electric fields 
was reduced below E\. The threshold electric field E\ 
decreases with increasing external magnetic field since 
the formation of immobile chains diminishes the effective 
spacing between cell plates resulting in higher electric 
field intensity for particles at the top parts of the chains. 
Upon exceeding E\ chains start to move and form lo- 
calized clusters of chains. At amplitudes of electric field 
above E2, the system transforms to a gas of short chains. 
Increase of the magnetic field above 50 Oe creates static 
long vertical chains touching the upper plate of the cell. 

We explored pure ac magnetic driving by placing 
the magnetic particles into a cell filled with toluene to 
dampen the kinetic energy in the system. Before each ex- 
periment, the system was driven to a gas-like state with 
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FIG. 3: Rescaled cluster area vs time at zero (a) and 10 Oe 
(b) magnetic fields. Solid lines represent simple functional 
dependence Q for different characteristic growth times r. 
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FIG. 4: Structures formed in external ac magnetic field: ring; 
compact clusters; chains of dipoles. Bottom: Patterns formed 
by nickel spheres (0.053 of the surface monolayer coverage) 
under magnetic driving at 20 Hz (clustered phase), 50 Hz 
(net-like structure), and 100 Hz (chains phase). 

an ac electric field in zero magnetic field to create a uni- 
form distribution of particles over the cell. Subsequently, 
the electric field was turned off and the system was sub- 
jected to ac 15 Oe magnetic field with frequency from 
to 200 Hz. The magnetic dipole-dipole interaction favors 
completely different particle organization. Fig. ^(upper 
part) demonstrates some selected structures. 

The nature of the self-assembled structures is quite 
straightforward. Since the local arrangement is domi- 
nated by a highly anisotropic dipole-dipole magnetic in- 
teraction, the head-to-tail interaction is the strongest, 
favoring formation of the quasi-one dimensional chain 
structures. Each chain has a certain rigidity against 
bending since in a bent chain the dipoles are not fully 
aligned. When the bending angle exceeds some thresh- 
old value, the energy barrier to bending starts to decrease 
since closing of the chain now becomes more favorable 
due to the attraction of the opposite ends of the chain. 
It can be shown that a ring is more energetically favor- 
able than a chain if the number of particles in the chain 
exceeds four 0. However, to form a ring from a chain, 
one needs to overcome a strong potential barrier associ- 
ated with chain bending. Consequently, a more plausible 
mechanism for forming a ring is a "fusion" of several 
short chain segments of appropriate orientation, polar- 
ization and position, see |S|- Chains can also form com- 
pact structure, such as a cluster of dipoles. Such dipolc 
clusters can nucleate in the situation when two neighbor- 
ing parallel chains have opposite polarization and attract 
each other to form a cluster. 

One can vary the frequency of the ac magnetic field 
to control the time it takes the particles (or chains) to 
rotate in order to keep its moment aligned with the local 
field and to move along the local field gradient. If the 
frequency is too high (200Hz in our case), nothing hap- 
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FIG. 5: Saturated chain length vs frequency of applied 15 Oe 
ac magnetic field for different amount of nickel 90/im particles 
in the cell. Solid lines are the fits to expression y — A+B/(x — 
fo). Inset: Saturated chain length vs applied ac (100Hz) 
electric voltage amplitude for magnetically driven (H max = 
15 Oe; f=25 Hz) system (0 ~0.10). 

pens since the characteristic reaction time of the system 
is higher than the period of the ac magnetic field. At 
some point while decreasing the frequency, the particles 
start to react to the external ac magnetic field and short 
chains appear in the system. A further decrease of the 
frequency leads to the creation of more energetically fa- 
vorable configurations such as rings and branched chains. 

Compact clusters as shown in Fig.0|can be initiated by 
decreasing the frequency of the ac magnetic field. The 
bottom panel of Fig. ^ shows various patterns formed 
under the influence of an external ac magnetic field at 
20, 50, and 100 Hz for a cell with 5.3% surface cover- 
age of Ni particles. Each experiment was started from 
randomly dispersed configuration of the nickel spheres 
over the bottom plate of the cell. Clearly, the pattern 
is determined by the frequency of the external ac mag- 
netic field: low-frequency excitations (0-30 Hz) produce 
a clustered phase; high-frequencies (80-200 Hz) assemble 
short chains while intermediate-frequencies favors net- 
like patterns. The resulting pattern is strongly history 
dependent. For example, a continuous decrease of the 
magnetic field frequency from 200 to 10 Hz results in the 
formation of a net-like structure, in contrast to the for- 
mation of clusters obtained when the system "relaxes" 
at a constant frequency of 10 Hz. Denser configurations 
(0=0.115 and 0.133), however, do not exhibit the clus- 
ter phase and instead, a phase transition to the network 
phase (infinitely long multi-branch chains) is observed at 
low frequencies. To demonstrate the phase transition to 
the network phase, the average chain length as a func- 
tion of the ac magnetic field frequency is plotted in Fig. 
[5] for different densities. Data was taken after the sys- 



tem was left to "relax" in the applied field for about 10 
minutes to attain their respective equilibrium state. As 
the frequency of the magnetic field decreases, the average 
length of the chains tends to diverge for dense configu- 
rations (0=0.115 and 0.133) and saturates for less dense 
ones. The solid lines in the figure are fits to the ex- 
pression, A + B/{x — fo), where f designates a critical 
frequency when the system undergoes a phase transition 
to the net-like phase. The critical frequencies extracted 
from the fits are 14.10 Hz, and 6.11 Hz for 0=0.133, and 
0=0.115, respectively. The critical frequency for surface 
coverage of 0=0.07 resulted in a negative value, / =-5.62 
Hz, indicating that there is no transition to the network 
phase. Indeed, for such low density configurations, com- 
pact clusters were observed at low frequencies. 

Inset to Fig. [3] demonstrates the effect of an ac electric 
field on the system of 240000 Ni particles (0 ~ 0.10) in ac 
magnetic field (H max =15 Oe; f=25 Hz). The frequency 
of the electric field was kept at 100 Hz to provide mostly 
two-dimensional motion. The average length of the chain 
decreases with increasing ac electric field, suggesting that 
the latter acts as an analogue to temperature. 

We studied the self-assembly of magnetic microparti- 
cles in ac electric and magnetic field. Excitation of the 
system by an ac magnetic field revealed a variety of pat- 
terns that can be controlled by adjusting the frequency 
and the amplitude of the field. We found that at low 
particle densities the low-frequency magnetic excitation 
favors cluster phase formation while high frequency ex- 
citation favors chains and net-like structures. For denser 
configurations an abrupt transition to the network phase 
was observed. We thank Ulrich Welp for help with mag- 
netization measurements. This research was supported 
by the US DOE, grant W-31-109-ENG-38. 
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